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In this paper we report numerically observed spontaneous vanishing of mean curvature on a 
' developable cone made by pushing a thin elastic sheet into a circular container 0. We show that 

this feature is independent of thickness of the sheet, the supporting radius and the amount of 
deflection. Several variants of developable cone are studied to examine the necessary conditions 
that lead to the vanishing of mean curvature. It is found that the presence of appropriate amount 
| of radial stress is necessary. The developable cone geometry somehow produces the right amount 

. of radial stress to induce just enough radial curvature to cancel the conical azimuthal curvature. 

Fj^ ■ In addition, the circular symmetry of supporting container edge plays an important role. With an 

elliptical supporting edge, the radial curvature overcompensates the azimuthal curvature near the 
minor axis and undercompensates near the major axis. Our numerical finding is verified by a crude 
experiment using a reflective plastic sheet. We expect this finding to have broad importance in 
describing the general geometrical properties of forced crumpling of thin sheets. 
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PACS numbers: 46.70.De, 68.55.Jk, 46.32.+X 

I. INTRODUCTION 



The crumpling of a thin sheet can be understood as the condensation of elastic energy into a network of two types of 
singular structures: folding ridges and point-like vertices. Scaling laws governing the energy and size of the ridge have 
Q ■ been obtained analytically and tested numerically [2-5]. Point-like singularities are also studied extensively [1,6-12]. 
Besides the understanding of these singularities, however, there is a lack of knowledge about the geometrical and 
mechanical effects of confining forces on a thin elastic sheet. 
' In this paper, we report an interesting observation of the response of an elastic sheet to the exertion of confining 
forces, using the developable cone or "d-cone" geometry studied by Cerda et al [8|[l([i3- Specifically, this system is 
£SJ ■ defined as follows. We push the center of a circular elastic sheet of radius R p axially into a cylindrical container of 
^sO ', radius R by a distance d. It is useful to express the deflection of the sheet by e = d/R. Due to the constraint of 
unstretchability, some part of the sheet is buckled inwardly and the sheet deforms into a non-axisymmetric conical 
surface that is only in partial contact with the edge of the container (FIG.^), thus forming a single developable cone. 
• In the limit of pure bending, that is, as thickness goes to zero, the shape of c£-cone has been shown to follow 
the solutions of the Elastica equation 13]. For asymptoticaly small deformation e, the shape can be analytically 
determined: this indicates that the buckled region occupies a fixed azimuthal angle of 138°, independent of all 
relevant length scales |l]|l3- However, for a physical sheet with finite thickness h, as discussed in Ref. [l^. stretching 
must occur and the stretching energy is concentrated in a small core region near thepushing tip. Outside this core 
1 ■ region, bending energy strongly predominates over stretching energy. Cerda et al |l| propose that the size of the 
core region scales as /1 1 / 3 ./? 2 / 3 . This scaling law was verified by our numerical simulations through direct geometrical 
measurements and through inferences from the measurements of central pushing force However, our analytical 
estimates indicate that this scaling cannot by asymptotic. Away from the core region, the confining normal forces 
from the edge of the cylindrical container arise from the requirement to counteract the central pushing force. In 
the situation of pure bending, Cerda and Mahadevan [l^ predict that the sheet exerts a concentrated point force 
to the container at the two take-off points where the sheet bends away from the container. Due to the translational 
symmetry along the region of the surface in contact with the container, the normal force is found to be independent 
of the azimuthal angle. 

The problem of interest here is the response of the sheet to these azimuthally constant normal forces from the 
edge of the container. Through our numerical study, we find that these normal forces produce a striking effect: 
they cause the mean curvature to vanish near the supporting edge within our numerical precision. In this work, we 
investigate this finding in a single d-cone as well as in variants of ci-done chosen to explore the conditions for this 
striking phenomenon. The numerical method we used is specified in Section II, consistent with that of Ref. ^3 • We 
also describe our validation of its accuracy. In Section III, we give the detailed description of our numerical findings 
and effects of altering the systems in various ways. In Section IV, a simple experiment is performed to verify our 
numerical observations for d-cone. We conclude with Section V and discuss limitations of and implications from our 
findings. 



o 



- 1—1 

X 



2 



A Z 




y 



FIG. 1: A simulated developable cone, formed as we axially push a hexagonal elastic sheet into a cylindrical container. The 
force F is applied to the central point of the sheet in the z direction. The dots are the lattice points of our simulation, as 
described in Section II. The buckled region appears at the lower right region of the sheet, where it bends downward and away 
from the container. 



II. NUMERICAL METHODS 

We begin by specifying the numerical models we use. An elastic sheet is modelled by a triangular lattice of springs of 
un-stretched length a and spring constant fc, after Seung and Nelson |14j . Bending rigidity is introduced by assigning 
an energy of J(l — fi\ - h-i) to every pair of adjacent triangles with normals h\ and fi2- When strains are small compared 
to unity and radii of curvature are large compared to the lattice spacing a, this model bends and stretches like an 
elastic sheet of thickness h = ay^8J/k made of isotropic, homogeneous material with bending modulus k, = J\/3/2, 
Young's modulus Y = 2ka/h\/3 and Poisson's ratio v = 1/3. Lattice spacing a is set to be 1. The shape of the sheet 
in our simulation is a regular hexagon of side length R p . The typical value of R p is 60a. 

To obtain a single d-cone shape, we need to simulate the constraining container edge and pushing force. As shown 
in FIG. m the edge lies in the x — y plane and is described by equation x 2 + y 2 = R 2 . Pushing in the center of the sheet 
is accomplished by introducing a repulsive potential of the form U{ mcc (zi) = —Fz\, where z\ is the z coordinate of the 
lattice point in the center and F is the magnitude of the pushing force. This force is applied in the p ositive z direction. 
The constraining edge is implemented by a potential of the form U e d ge = C p H(zi)/([(^/x 2 + y 2 — R) 2 + zf ] 4 + £ 8 ), 
where £, C p are constants and the summation is over all lattice points with coordinates (xi,yi, Zj). H{z) is the unit 
step function smoothed over a lattice spacing, which makes certain that this potential only acts on the lattice points 
that have already moved into the container (those with Zi > 0). We choose the range £ of the potential to be one 
lattice spacing. Thus the force decays rapidly once the lattice points go away from the edge. A discussion concerning 
the choices of the edge potential form is carried out in Section V. 

The conjugate gradient algorithm |17| is used to minimize the total elastic and potential energy of the system as a 
function of the coordinates of all lattice points. This lattice model behaves like a continuum material provided that 
the curvatures are everywhere much smaller than 1/a. This limitation restricts the values of deflection e of d-cone to 
be below 0.25. 

We determine the curvatures approximately from each triangle in the sheet. For this measurement, we take the 
curvature tensor to be constant across each triangle. We calculate it using the relative heights of the six vertices of 
the three triangles that share sides with the given triangle |5j. The six relative heights Wi normal to the triangle 
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surface are fit to a function of the form 

w l = b\ + b 2 u t , + b^Vi + b^ul + b 5 UiVi + b 6 vf , i = 1, . . . , 6 (1) 

where {ui,Vi,Wi} are coordinates of the vertices in a local coordinate system that has w axis perpendicular to the 
surface of the given triangle. This choice of local coordinate system ensures that 62 and 63 are negligible so that 
curvature tensors can be determined only from the coefficients of quadratic terms. In practice, our numerical findings 
do show that the values of 62 and 63 are on the order of 1CP 2 or lower. Therefore, curvature tensors follow immediately 
from the identification C uu — 2 x b^.C vv — 2 x b$,C uv = 65. Mean curvature C is defined as half of the trace of 
curvature tensor: C = (C uu + C vv )/2. 

This lattice model of elastic sheet has been used to study both the ridge and the point-like singularities in crumpled 
sheets [2-5] [H^j- The accuracy of it has been tested in various ways. Using this model, Lobkovsky Q numerically 
verified the "virial theorem" for ridges that bending energy is five times the stretching energy, for asymptotically 
thin sheets. In the study of d-cone, the ratio of the normal force from the (^-function term to that from the angle- 
independent term is found to be 0.70 numerically |12f. compared well with the theoretical prediction 0.69 Also, 
the azimuthal profiles of curvature are in reasonably good agreement with the theoretical prediction as deformation of 
d-cone goes to zero ^| . In addition, we tested our program by setting the lattice at different initial states and letting 
the program look for minimized energy state. We found the lattice converges to the same state, with agreements of 
pushing force and energies better than one percent for fixed deflection. 



III. CONDITIONS FOR VANISHING MEAN CURVATURE 



We consider the mean curvature of the sheet along a radial line that lies in the non-buckled region. Ideally, for a 
regular conical surface, the mean curvature only comes from the azimuthal component Cgg and follows a 1/r decay, 
where r is the distance to the tip. For a <i-cone described here, however, due to the pushing of the normal forces 
from the edge of container, the sheet has to experience a small deformation near the region where it touches the edge. 
More specifically, these normal forces cause a small inward deflection of the sheet and hence induce both curvature 
and strain in the contacting region. Noticing that the curvature induced is in radial direction and has opposite sign 
to Cgg, we expect the mean curvature to be reduced near the edge. 
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FIG. 2: Radial profiles of mean curvature in three different directions on a d-cone surface. The d-cone has thickness h = 0.102a, 
confining radius R = 38a and deformation e = 0.15, where a is the lattice spacing in the numerical model [l4|. Inset shows the 
shape of the d-cone, where the three radial directions are denoted by the same symbols as the corresponding data. The slopes 
of fitted lines are —0.934, —0.936 and —0.820 for data denoted by pluses, crosses and circles, respectively. 
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Numerically, FIG. shows the mean curvature profiles along three different radial directions on the same (i-cone 
with thickness h = 0.102a, deformation e = 0.15 and cylinder radius R — 38a. Data denoted by pluses and crosses are 
measured in the non-buckled region of the sheet; data denoted by circles are measured in the buckled region. Inset 
shows the d-cone shape and the three radial directions, denoted by the same symbols as the corresponding data. The 
values of three fitted slopes are all close to —1, verifying the 1/r feature of the mean curvature. However, for profiles 
in the non-buckled region, an abrupt drop of the mean curvature is observed near r = 38a, where the sheet touches 
the edge of container. The curvature drops by more than one order of magnitude, making the mean curvature at 
the valley effectively zero. This indicates that the radial curvature caused by the pushing of normal forces from the 
edge not only reduces the mean curvature near the edge, but also virtually equals the original azimuthal curvature in 
magnitude. 

We further observe that the feature of vanishing mean curvature does not depend on particular values of R or e. 
FIG.Elshows the radial profiles of mean curvature for d-cone formed with four different confining radius R, with fixed 
thickness h = 0.102a and deformation e = 0.15. We observe that mean curvature drops down to zero effectively at 
each of the four radii. Moreover, FIG. 01 displays the mean curvature profiles along a radial line on d-cone formed 
with four different deformations, with fixed thickness h = 0.102a and confining radius R — 38a. It is seen that for 
four deformations, the mean curvatures all go to zero near r = 38a, which the sheet touches the edge. 
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FIG. 3: Radial profiles of mean curvature for d-cones with four different values of confining radius: R = 26a, 32a, 38a, 44a. 
Thickness h — 0.102a and deformation e = 0.15. 



To explore the conditions for the vanishing of moan curvature, we study variants of d-cone and make comparisons. 
To determine what is important in this phenomenon, we look at the the balance of normal force on each element of 
the sheet, which is expressed by the force von Karman equation |16| 

d a dpM af3 = a a pC af3 + P . (2) 

here M a p are torques per unit length, a a p are in-plane stresses, and P is the external normal pressure on the sheet. 
For an asymptotically thin sheet, the azimuthal curvature of d-cone is Cgg = <j)(0)/r, where (f)(9) — — e for \9\ > 9 
(unbuckled region). In the limit of small deformation, it is found that (f)(9) = —ecos(a9)/ cos(a9o) for \9\ < 9q (buckled 
region), with 9q — 1.21 rad and a = 3.8 The torques are proportional to curvatures through a constitutive 

law implicit in the energy equations ^^:M rr = kvCqq — nv(f)/r and M$$ — kCqq = K(f)/r, where k is the bending 
modulus and v is the Poisson ratio. Using these laws to eliminate torques from the force von Karman equation, one 
obtains 



&90 = 




where the dots above functions denote 9 derivatives. For the unbuckled region, cf>/<j> = 0, so age = 2nv/r 2 . For 
the buckled region, <p/(f) — —a 2 , so age = —n(a 2 — 2v)/r 2 . Therefore, the azimuthal stress is compressive in the 



5 



10" 



O 



° ° 8 S g g x „ x x x 

I I I + + + ° °°°!°° 5s 3gSS66gS3o>< 



xxxxxx^ 



r x e=0.20 
o 8=0.15 

+ £=0.10 

* £=0.05 



30 38 



FIG. 4: Radial profiles of mean curvature for d-cones with four different values of deformations: e = 0.05, 0.10, 0.15, 0.20. 
Thickness h = 0.102a and confining radius R = 38a. 



buckled region but tensile in the unbuckled region. Cerda and Mahadevan show that for the same shape formed 
under a conical support, unlike a d-cone formed under a ring support as illustrated above, the azimuthal stress in the 
unbuckled region is compressive due to the exertion of normal pressure everywhere. In our case, the tensile azimuthal 
stress in the unbuckled region is changed to be compressive only near the edge of container, where the normal forces 
are nonzero. Whether the stresses are tensile or compressive, their magnitude goes as l/r 2 . We call them type I 
stresses and denote them by symbol . It is noted that the type I stresses exist both in azimuthal and radial 
directions. 

Having looked at the local normal force balance of the surface, we now investigate the global force balance. We 
consider a region that encloses area between inner radius R c and outer radius r, where r can take values between 
R c and R. R c denotes the size of the stretching core region. The tension exerted on the inner edge of this region is 
equivalent to the central pushing force of a regular <i-cone. This force must be balanced by the force due to radial 
in-plane stress a rr on the outer perimeter of the region. Let (3 be the angle between a radial line or generator of the 
d-cone and the horizontal. Since tan j3 = (f)(0), we have sin /3 — (f)(0) / y/l + (f) 2 (0). The balance of vertical forces yields 



J a rr r sin (3d0 — F , (4) 

which holds for every value of r from R c to R. It is easy to see that the type I stresses alone can't satisfy this equation, 
since they go as l/r 2 , they would give a l/r prefactor on the left side of equation, while the right-side of equation 
is independent of r. Therefore, the integral from type I stresses must vanish and there must exist some additional 
tensile stresses in the outer region that scale as l/r to satisfy Eq. Q}. We call these stresses as type II stresses 
and denote them by . They persist up to the supporting container edge, where normal force from the container 
counteracts the external pushing force. We can write er( 2 V = Fe(0), where e(0) is a function only of and satisfies 
J e(0)sm (3d0 — 1. It is obvious that e(0) is of order unity. To estimate the magnitude of type II stresses, we notice 
that F — dE/dd = (dE/de)/ R « k/R. Thus « F/r w n/(rR). The type II stresses are comparable with type I 
stresses only near the container edge since the ratio a^/a^ w r/R for R c < r < R. 

Now that we see that the contacting normal force and the stresses seem to matter, we vary the system in ways that 
alter them. For this purpose, we study several geometries that we call the cut cone, regular cone, puckered cylinder 
and elliptical d-cone, and compare them with regular d-cone. 

First, to study the effect of stresses, we cut the sheet along a radial line and push it into the same container to form 
a "cut cone" . Now that the boundaries along the cut line are free to move, the buckled region no longer appears and 
the sheet will necessarily overlap itself in order to fit into the container, as illustrated in FIG. [5j The cut cone mostly 
follows a conical shape as the unbuckled region of d-cone does, indicating that the azimuthal tensile stress remains 
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the same as d-cone. Near the free boudaries, however, the sheet deviates from a conical shape and does not contact 
against the confining edge. Instead, there is a section of the sheet near the free boundaries that has less curvature 
than the corresponding d-cone shape, due to the lack of bending moment there. Consequently, the free boundaries 
contact the container at a nonzero angle and the azimuthal stresses arise from the pushing force from the edge. Due 
to the decreases in the energies, the pushing force is reduced by a finite factor relative to that of d-cone. This causes 
the radial tensile stress to be reduced by the same factor. Numerically, it is found that the pushing force is reduced 
by about one order of magnitude. The measurements of the mean curvature along radial lines show rather small 
decreases in the contacting region, far from enough to make mean curvature vanish. These contrasts to the d-cone 
show that the appropriate amount of radial stress in the d-cone may play an important role for the mean curvature 
to vanish. 




FIG. 5: After we cut the sheet along a radial line, the buckled region no longer appears as we push the sheet into the same 
container. Instead, the free boundaries overlap to accommodate the confining edge. 



Next, we make a regular cone by removing one sixty-degree sector of the hexagon and joining the two free sides 
(see FIG. 0- We then push this regular cone into the same container as before, and measure the corresponding radial 
profile of mean curvature. The cone does not buckle for the range of forces we consider. We find that as we push 
with the similar force as that required for a d-cone of the same thickness and deformation, we observe the similar 
feature of vanishing mean curvature near r = R. To illustrate this point, FIG. [7| shows the radial profiles of mean 
curvature of a regular cone pushed with three different forces against the confining edge. The pluses, circles and 
triangles are for a regular cone pushed by half of the corresponding d-cone force, one (i-cone force and twice the 
d-cone force, respectively. We observe that the mean curvature vanishes as we push with one d-cone force, while it 
is not vanishing for half (i-cone force and is changing sign for double d-cone force. When the central pushing force is 
similar, the normal forces from the edge are similar for regular cone and d-cone. So are the azimuthal compressive 
stresses induced by them in the contacting region. These facts show that the pushing force required to make d-cone is 
somehow just enough to produce the right amount of stresses for which the induced radial curvature cancels original 
azimuthal curvature. 

We next investigate a variant that preserves the feature of a buckled sheet caused by a constraining boundary 
without the d-cone's feature of an additional force applied away from this boundary. To see this, we consider a 
puckered cylinder formed by rolling a sheet of cylindrical shape with radius R into a circular cylinder of radius b < R, 
as studied by Cerda and Mahadevan [l^. To make closer comparison with the d-cone, we replace the cylindrical 
confining surface with a circular ring of the same radius. Now that the puckered cylindrical sheet is confined by the 
ring and the ring exerts normal forces on the sheet, we expect the mean curvature to be reduced in the contacting 
region where the sheet touches the confining ring. The simulated shape is shown in FIG. 00 We measure the mean 
curvature along a longitudinal line parallel to the axis of cylinder. A small drop of mean curvature is observed where 
the sheet touches the ring. Not surprisingly, the reduction of mean curvature goes up as the ratio R/b increases. 
However, we find that the drop is no more than 10% for various values of R and b that can produce reasonable 



FIG. 6: A simulated regular cone formed by removing one sector and sewing the free boundaries. 
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FIG. 7: Radial profiles of mean curvature of a regular cone pushed with three different forces against the confining edge. The 
regular cone has thickness h — 0.127a and deformation e = 0.66. Confining radius is R — 38a. Plus is for regular cone pushed 
with half d-cone force; circle is for regular cone pushed with one d-cone force; triangle is for regular cone pushed with double 
d-cone force, with downward triangles denoting mean curvatures changing sign. 



confined shape, even for R/b = 2. The reason that the mean curvature doesn't vanish near the contacting region is 
attributed to the absence of stess in the non-azimuthal direction. By contrast, the d-cone has a radial tensile stress 
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which transmits the central force to the container. 




FIG. 8: A simulated puckered cylinder confined by a circular ring. 



Lastly, we observe that the circular symmetry is important for cancellation of mean curvature. To see this, an 
elliptical confining edge is used instead of a circular edge. Our original hexagonal sheet is pushed against this 
elliptical edge. The lengths of the major and minor axes of the ellipse are 45a and 30a, respectively, with a being the 
lattice spacing in the numerical model. FIG. |S] shows mean curvature profiles in two different radial directions, both 
of which lie within the unbuckled region of the sheet. One profile is measured along a major axis of the elliptical edge. 
The other profile is measured along the lattice direction sixty degrees from the first and thus meets the ellipse close to 
the minor axis. For the sake of uniform comparison, we measured curvatures only along principal lattice directions. 
It is seen that the first profile shows a drop of mean curvature near the contacting region, although this drop is not 
as pronounced as those in FIG. [3 On the other hand, the second profile shows that the mean curvature changes 
sign. The induced radial curvature is greater than the azimuthal curvature in magnitude, thus overcompensating the 
azimuthal curvature. We conclude that the break of circular symmetry in the unbuckled region leads to nonuniform 
degrees of compensations of mean curvature in different directions. 



IV. AN EXPERIMENTAL VERIFICATION 



To demonstrate that our numerically observed phenomenon occurs also in physical sheets, a crude experiment was 
performed. First, we made a d-cone by pushing a reflective plastic sheet into a circular container of radius 9.1cm. 
The sheet is a an ordinary 22cm x 28cm projection transparency used in presentations. Its thickness is about 0.1mm. 
We measured the curvature profile of this sheet by observing the images of standard objects reflected in the sheet. 
The magnification of each object allows us to infer the curvature at the corresponding position on the d-cone. The 
experimental setup is sketched in FIG. EH The standard objects are identical black circles arranged in a row, printed 
on a second transparency. The row of circles is oriented such that its images on the surface of the d-cone lie along 
a radial line in the unbuckled region. The images of the circles are turned into ellipses that are progressively wider 
in the azimuthal direction as they approach the vertex, since the azimuthal curvature is inversely proportional to 
the distance to the vertex. In the noncontacting region, the images all have the same height in the radial direction, 
since radial curvature is zero there. However, the image on the edge is different. It is being contracted in the radial 
direction due to the induced radial curvature there. If the radial curvature is of the same magnitude as azimuthal 
curvature at the edge, the contraction in the radial direction and the expansion in the azimuthal direction have to 
be by the same factor. Let w be the azimuthal size and I be the radial size of the image, then the w/l ratio for the 
image at the edge has to be the square of that ratio for the image just outside or inside the edge. 
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FIG. 9: Radial profiles of mean curvature in two different directions for the same d-cone confined by a elliptical edge with 
major axis of 45a and minor axis of 30a. The distance is in linear scale. The data denoted by circles are measured along a 
major axis of the elliptical edge (long direction). The data denoted by pluses and crosses are measured in the direction that 
makes 60 degrees to the major axis (short direction). The crosses denote that the mean curvature becomes negative, so that it 
is overcompensated. 




FIG. 10: (Color online) Sketch of apparatus for measuring curvature. Reflective transparency A is pressed into the container 
B by a probe C. Object D, a row of identical black circles on another transparency, forms the reflected image E, which is 
photographed by camera F. 



We used a digital camera to take pictures of the images near the edge. FIG. ^2 is one of the pictures taken for 
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TABLE I: The sizes of images in FIG. 1111 from top to bottom, in arbitrary but fixed units, w is the azimuthal size (width) 
of the image; I is the radial size (height). The third spot lies on the edge. The actual size of the top spot is about 1.5mm in 
width and 1mm in height. 
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TABLE II: The sizes of images in FIC I12I from top to bottom, in arbitrary but fixed units, w is the azimuthal size (width) of 
the image; I is the radial size (height). The top one lies on the edge. The actual size of the top spot is about 3mm in width 
and 2mm in height. 

e w 0.25. It shows images of 5 circles that lie along a radial direction. The lower images are closer to the tip. The 
central spot is the one that lies on the edge. The values of azimuthal and radial sizes of the images as well as the their 
ratios are listed in table I, with the corresponding ratios that show comparison highlighted. From the table, we can 
see that the expected relationship between w/l ratio of the image that lies on the edge and that of the adjacent image 
is favored by the data, indicating that radial curvature is indeed of the same magnitude as azimuthal curvature. As 
another check, FIG. El is the picture taken for e w 0.15. It shows images of 7 circles. The images are a little out of 
focus, and we measure their sizes according to the black core region, excluding the blurry fringes. The top spot lies 
on the edge. Table II gives the values of their sizes and corresponding ratios. We observe that these data are also in 
agreement with our numerical observation. 



V. DISCUSSION AND CONCLUSION 

In this paper, we have presented our findings of vanishing mean curvature in the contacting region in a single 
developable cone. The mean curvature near the supporting edge is observed to vanish within our numerical precision 
for d-cones with various confining radii, thicknesses and deformations. We investigate this phenomenon by studying 
variants of d-cone and identifying the features that are important for the vanishing of mean curvature to happen. 
It is found that the presence of appropriate radial stress is necessary. We also find that the d-cone in a circular 
container somehow produces the right amount of normal forces that induce just enough radial curvature to cancel 
azimuthal curvature. By comparing <i-cones formed with circular and elliptical confining edges, we find that circular 
symmetry is indispensable for the mean curvature to vanish. The numerically observed phenomenon is verified by a 
crude experiment of projecting identical circles onto the d-cone surface and comparing the corresponding magnification 
ratios. 

It is worth noting that the cancellation of curvatures does not depend on particular choices of edge potential. In 
our assessment as long as the choices of potential represent a sharp edge, as is the case for our choices presented in 
Section II, the induced radial curvature cancels the azimuthal curvature at the edge. Some supporting evidence is as 
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FIG. 11: (Color online) Images of circles as reflected from the surface of a d-cone with e m 0.25. The lower images are closer 
to the tip. The third image lies on the edge, which is perpendicular to the line of images. 



follows. If the radial curvature were induced due to the particular form of the edge potential, that radial curvature 
would be independent of the deflection of the sheet. However, this contradicts with our observations that induced 
radial curvature does depend on the deflection of the sheet. 

The vanishing of mean curvature is precisely the property of a "minimal surface", i.e. a surface of minimal surface 
area for given boundary conditions. Liquidlike membranes, such as a soap film |l8| . follow the shape of a minimal 
surface in equilibrium and thus have zero mean curvature everywhere. This phenomenon is purely geometric and 
it reminds one of the Gauss-Bonnet theorem, which constrains the average Gaussian curvature in a surface jlflj . 
Our results for elliptical container suggest that it is some average mean curvature along the rim that vanishes. 
Understanding radial curvature C rr requires more than Elastica equation, since C rr necessarily creates stretching 
not allowed in the Elastica approach |13| . One approach to understanding the C rr is to assume the contact force 
inferred from the Elastica approach and use the force von Karman equation to infer C rr . This unexpected geometrical 
regularity in force thin sheets is important to understand, since it surely influences a broad range of situations extending 
far beyond those studied here. We shall explore this phenomenon further in future publications. 
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FIG. 12: (Color online) Images of circles as reflected from the surface of a d-cone with e m 0.15. The lower images are closer 
to the tip. The first image from the top lies on the edge, which is perpendicular to the line of images. The images are a little 
out of focus. The sizes are measured based on the black core region, excluding the blurry fringes. 
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